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ORIGINAL ARTICLE 

Nonlinear Population Pharmacokinetics of Sirolimus 
in Patients With Advanced Cancer 

K Wu^ EEW Cohen^ 3^ LK House\ J Raml'rez^ W Zhang^ MJ Ratain^ 2,3 rr 51^34 

Sirolimus, the prototypical inhibitor of the mammalian target of rapamycin, has substantial antitumor activity. In this study, 
sirolimus showed nonlinear pharmacokinetic characteristics over a wide dose range (from 1 to 60 mg/week). The objective 
of this study was to develop a population pharmacokinetic (PopPK) model to describe the nonlinearity of sirolimus. 
Whole blood concentration data, obtained from four phase I clinical trials, were analyzed using a nonlinear mixed-effects 
modeling (NONMEM) approach. The influence of potential covariates was evaluated. Model robustness was assessed using 
nonparametric bootstrap and visual predictive check approaches. The data were well described by a two-compartment 
model incorporating a saturable Michaelis-Menten kinetic absorption process. A covariate analysis identified hematocrit 
as influencing the oral clearance of sirolimus. The visual predictive check indicated that the final pharmacokinetic model 
adequately predicted observed concentrations. The pharmacokinetics of sirolimus, based on whole blood concentrations, 
appears to be nonlinear due to saturable absorption. 
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Sirolinnus (rapamycin) is commonly prescribed as an immu- 
nosuppressant and is indicated for the prevention of allograft 
rejection J As the prototypical inhibitor of the mammalian tar- 
get of rapamycin, sirolimus (like other mammalian target of 
rapamycin inhibitors) has substantial antitumor activity both 
in animals and humans.^-^ Sirolimus has low bioavailability 
(14% on average) and a long terminal elimination half-life 
of -62 h/ Studies in transplant patients have demonstrated 
marked interindividual pharmacokinetic variability, resulting 
in the widespread use of therapeutic drug monitoring based 
on whole blood concentrations. ^-^^ 

There are several publications regarding the pharmacoki- 
netics of sirolimus in transplant patients, although none have 
explicitly incorporated the nonlinear pharmacokinetic charac- 
teristics into a mixed-effects population model.^"^^ Jiao et al.^^ 
reported the nonlinearity in the pharmacokinetics of sirolimus, 
but they were unable to develop an explicit model to describe 
this due to limited measurements. Recently, several phase I tri- 
als of sirolimus in patients with cancer were completed, includ- 
ing the intermittent administration of higher doses.^"^'^^ These 
studies provided the opportunity to investigate the nonlinearity 
in sirolimus disposition using whole blood concentration mea- 
surements. The detection and characterization of nonlinearities 
provided by population modeling allows a better understand- 
ing of how a drug should be used in clinical practice.^^ 

The objective of this study was to develop a population 
pharmacokinetic (PopPK) model for sirolimus, while exploring 
possible nonlinear absorption characteristics in whole blood 
measurements. These measurements were obtained from 
clinical trials of patients with advanced cancer who received 
sirolimus in a wide range of dosages (1-60 mg/week). This is 
the first PopPK report of sirolimus in patients with advanced 
cancer. 



RESULTS 

A total of 563 concentration data points from 76 patients 
with advanced solid tumors enrolled in four different phase 
I trials at The University of Chicago were available for the 
analysis.^^'^^ An example of the data records is provided in 
the Supplementary Table S1 online. 

Noncompartmental analysis 

Figure 1a,b demonstrates the nonlinearity of sirolimus phar- 
macokinetics. The slope of dose-normalized area under 
the curve (AUC)q_ vs. dose differed significantly from zero 
(P< 0.01), indicating that the drug exposure did not increase 
linearly with dose. The terminal elimination half-life did not 
deviate significantly from zero (P > 0.05), suggesting that the 
elimination is linear over the dose range in this study. 

Nonlinear PopPK model 

Sirolimus concentration vs. time curves were best described 
using a two-compartment model with a saturable absorption 
model (Michaelis-Menten equation) (Figure 2). The base 
model was characterized by the following expressions: 
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where A^, A^, and A^ are the amounts of drug in the intestinal 
lumen and central and peripheral compartments. Vm (|jg/l-h) 
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Figure 1 The ratio of area under the curve (AUC)q_ and dose 
(AUCQ_/dose) (a) and terminal elimination half-life (b) vs. doses. 
Males and females are represented by dots and squares. The data 
from trials 2 (dose range 1-16 mg) and 3 (dose range 15-35 mg) 
are shown in pink and blue. 



where C^^^ and C^^^^ are the observed and predicted concen- 
trations. £^ and £^ are randonnly distributed variables with a 
nnean of zero and variances of and accounted for the 
residual variabilities. 

The estinnated PopPK paranneters are listed in Table 1 .The 
relative standard errors of paranneter estinnates ranged fronn 
8.1 7% to 55.4%. Figure 3 shows the relationship between the 
observed and population nnodel-predicted concentrations and 
the relationship between the observed and individual nnodel- 
predicted concentrations. The subjects in these studies were 
outpatients, and undocunnented variability in tinning of doses 
nnay be a significant factor innpacting both bias and preci- 
sion.^^ The MS nnethod was tried as there were 16 sannples 
(2.0% of all observations) below the linnit of quantification.^^ 
It was not included in our final nnodel because addition of the 
M3 nnethod did not innprove the nnodel fit. Hennatocrit was the 
only significant covariate affecting the apparent clearance of 
sirolinnus (clearance decreased with increasing hennatocrit). 
Drug fornnulation (liquid vs. solid) did not have a significant 
innpact on the absorption-related paranneters. 

Model evaluation 

The nnedian paranneter values resulting fronn the bootstrap 
procedure agreed with the estinnates fronn the final popula- 
tion nnodel. This suggests that the paranneters in the final 
nnodel were reasonably well deternnined and the nnodel 
was stable. Fronn 1 ,000 bootstrap runs, 985 nnininnized suc- 
cessfully and were included in the bootstrap analysis. The 
results of the bootstrap analysis are sunnnnarized in Table 1. 
Figure 4 shows the nnedian and the 5th and 95th percent 
prediction intervals fronn the visual predictive check sinnula- 
tion with the observed data superinnposed. These plots show 
that nnost of the observed concentrations on all dose levels 
fell within the 5th-95th percent prediction interval. Observed 
concentrations <10% lay outside the prediction intervals. The 
visual predictive check shows that the final nnodel adequately 
describes the nnajority of the data. 
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Figure 2 Two-compartment model with a saturable Michaelis- 
Menten kinetic absorption. CL^, clearance from central compartment; 
CL^, clearance between central and peripheral compartments; Gl, 
gastrointestinal tract; Km, the amount at 50% of Vm; Vm, maximum 
absorption rate. 

is the nnaxinnunn absorption rate and Km (nng) is the drug 
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and V^ were 



fitted. At 0 h, A^ 

A connbined proportional and additive error nnodel was 
used to describe the residual unexplained variability repre- 
senting the variance between the observed concentrations 
and those predicted by the nnodel: 



^obs=^pred-(1 +^i)+^2 



DISCUSSION 

This is the first PopPK study of sirolinnus in patients with 
cancer. Using whole blood concentrations, the data were 
adequately explained by a Michaelis-Menten absorptive pro- 
cess, which results in increased apparent oral clearance with 
dose. 

In addition, our nnodel dennonstrated that the appar- 
ent oral clearance of sirolinnus was inversely associated 
with hennatocrit, although this has only a nnodest effect on 
clearance (e.g., an increase in hennatocrit fronn 35.1% to 
44.7% results in a decrease in apparent clearance fronn 
12.9 to 12.4 l/hr). Sirolinnus is distributed nnostly into red 
blood cells and very little into plasnna (<5%).^^-^^ The 
sequestration of sirolinnus, like tacrolinnus, into blood cells 
is believed to be nnainly due to the presence of binding pro- 
teins such as FK binding protein. Lipid solubility, degree 
of ionization, nnolecular size, and hydrogen-bonding ability 
have been identified as the nnain deternninants for uptake 
by the red blood cells. The relationship between clear- 
ance and hennatocrit nnay reflect the exchange of sirolinnus 
between plasnna and red blood cells. This exchange plays 
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Figure 3 Observed (DV) vs. (a) population model predicted and (b) individual model predicted concentrations. The solid lines are diagonal 
lines of identity, (c) Conditional weighted residuals vs. population predictions. 



Table 1 Final population pharmacokinetic parameter estimates and bootstrap 
results of sirolimus 



NONMEM 



Bootstrap 



Parameter 

(unit) Definition 



Estimate (% SE) 
(shrinkage) 



Estimate 
(% SE) 



CL/F= 6^ X (median/hematocrit)^ 
0, (l/h) 



0. 

V,/F{\) 
CL,/F(l/h) 

VJF{\) 

Vm (|jg/l-h) 
Km (mg) 



Typical value of 
clearance from central 
compartment 

Factor of hematocrit 

Volume of distribution 
of central compartment 

Clearance between 
compartments 

Volume of distribution 
of peripheral 
compartment 

Maximum absorption 
rate 

The amount at 50% 
of Vm 



12.9 (16.3%) 

0.14(55.4%) 
53.4 (38.0%) 

29.0 (8.17%) 

611 (11.3%) 

4.56 (37.7%) 
13.8 (50.3%) 



Interindividual variability (%SE) (shrinkage) 
CL/F 52.4% (57.8%) (1 2.5%) 

V^/F 52.4 % (57.8%) (24.0%) 

CLJF 70.5% (17.2%) (32.4%) 

VJF 1 9.3% (1 77%) (28.0%) 

Intraindividual variability (%SE) (shrinkage) 
Proportional (%) 2.17% (97.0%) (11.1%) 

Additive (ng/ml) 0.5 (35.5%) (1 1 .9%) 



12.8 (49.9%) 

0.11 (22.8%) 

59.1 (52.3%) 

27.8(10.1%) 
607(16.4%) 

4.61 (41.2%) 

14.2 (59.3%) 

66.7% (70.2%) 
66.7% (70.2%) 
69.3% (40.9%) 
20.5% (150%) 

2.53% (15.6%) 
0.48 (30.9%) 



The median of hematocrit is 35.1%. 
NONMEM, nonlinear mixed-effects modeling. 



an innportant role in its pharnnacokinetic behavior.^ Many 
other drugs, such as tacrolinnus^^ and cyclosporine,^^ also 
show a correlation between whole blood clearance and 
hennatocrit.^^'^^ 

In the forward addition process, gender was found to be 
significantly related with V^/F (\//Ffor fennale is 75% of 
the \//Ffor nnale), but did not pass the backward elinnina- 
tion. The volunne of distribution for the peripheral connpart- 
nnent increased slightly with increasing body weight (the 
exponent is 0.676); however, the addition of body weight 
to volunne of distribution as a covariate did not result in 
a statistically significant change in the objective function 
value (1.589). 

Michael is-Menten kinetics was incorporated in our nnodel 
to describe the saturable absorption. This is consistent with 
a previous study on intestinal absorption, which dennon- 
strated a concentration-dependent and saturable transepi- 
thelial transport of sirolimus across Caco-2 monolayers.^^ 
The mechanism of absorption is not very clear, and the 
model could be further understood once the drug transport 
mechanisms are elucidated. Cummins et al.^^ studied the 
cellular pharmacokinetics of sirolimus in CYP3A4-trans- 
fected Caco-2 cells. They found that the net effect of P-gly- 
coprotein on metabolism was not as great as that expected 
for sirolimus, a P-glycoprotein substrate, and assumed that 
there might be multiple transporters involved in sirolimus 
absorption. 

Our model fit the concentration data well. The residual 
of the pharmacokinetic parameter estimates was accept- 
able (Table 1). The clearance of central compartment in our 
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Figure 4 Visual predictive check of the model for all dose levels. The median model predictions are shown by solid lines and 95% prediction 
intervals are encompassed by the broken lines. The dots present the observed concentrations. 



study is connparable with those reported previously (7.9-28.3 
l/l^ J 10,12,20,23 apparent volunne of distribution of the periph- 
eral connpartnnent (tissue) is in the range of values reported in 
previous studies (1 1 .6-1 ,350 j). 8-10,12,20,23 Because sirolinnus is 
hydrophobic, it is widely distributed in the lipid nnennbranes of 



body tissues and the erythrocytes and shows a large appar- 
ent volunne of distribution. ^°'22 

The current PopPK analysis can be used to predict the 
effect of dose and/or schedule change on the whole blood 
concentrations of sirolinnus. This will facilitate future studies 



CRT: Pharmacometrics & Systems Pharmacology 



Nonlinear PopPK of Sirolimus in Patients With Advanced Cancer 

Wu etal 



of this agent in addition to providing information tliat may 
guide dosing for patients individually. 

METHODS 

Patients and data collection. The clinical trials were reviewed 
and approved by The University of Chicago Institutional Review 
Board. ^^'^^ Written informed consent was obtained from all 
patients. The demographic characteristics of the patients are 
summarized in Table 2. Treatment was administered on an 
outpatient basis according to the treatment schedule listed in 
Table 3. Patients self-administered sirolimus (without food) at 
the scheduled times per protocol. Patients from trial 1 were 
given high dose sirolimus. In trials 2 and 3, escalated doses of 
sirolimus were administered. In trial 4, a fixed sirolimus dose 
of 4 mg was given. Dosing was independent of body weight. 
Sirolimus was supplied in liquid form for trials 1 , 3, and 4, and 
in tablet form for trial 2. 

Blood sampling and sirolimus analysis. Sirolimus whole blood 
concentrations were obtained according to the sampling 
schedule shown in Table 3. Samples were stored at -80 °C 
until analysis by high-performance liquid chromatography com- 
bined with tandem mass spectrometry (LC/MS/MS) for trials 2, 
3, and 4. The extraction and LC/MS/MS methods are described 
in Supplementary Appendix online. The samples from trial 1 
were analyzed using high-performance liquid chromatography 
and LC/MS/MS. The limit of quantification (in ng/ml) was 2 
(trial 1), 0.28 (trials 2 and 3), and 0.49 (trial 4). 

Noncompartmental analysis. To evaluate the dose depen- 
dence of the pharmacokinetic behavior of sirolimus, AUCq_ 
and half-life were calculated using a noncompartmental 
pharmacokinetic analysis implemented in the PK Solutions 
software (version 2.0; Summit Research Services, Mon- 
trose, CO). Sirolimus whole blood concentrations from trials 
2 and 3 were used for this analysis. To establish dose linear- 
ity and proportionality, dose-normalized AUCq_ and half-life 
were analyzed using regression analysis (a = 0.05). 

Nonlinear PopPK model. An analysis was performed using a 
nonlinear mixed-effects modeling (NONMEM) approach as 
implemented in NONMEM (version VII, level 1; ICON, Ellicott 
City, MD) in conjunction with a gfortran compiler.^^ First-order 
conditional estimation with interaction and the ADVAN 13 sub- 
routine were applied. 

A number of candidate models were assessed in describ- 
ing the concentration data and apparent nonlinearity of the 
pharmacokinetics of sirolimus: two-compartment model with 
Michaelis-Menten elimination; Michaelis-Menten absorp- 
tion; zero-order absorption and first-order loss in gastrointes- 
tinal tract; parallel zero-order and first-order absorption; and 
Weibull absorption.^ Mixture distributions were also explored 
to examine whether the absorption patterns had a multimodal 
characteristic. Demographic data of the patients, shown in 
Table 2, were tested in the model one by one as candidate 
covariates on each of the parameters. The final PopPK model 
was established using the stepwise fonA/ard addition and back- 
ward elimination method.^^ Comparison of the models was 
based on the objective function value provided by NONMEM 
at a significance level of 0.05 (equal to a decrease of 3.84 in 



Table 2 Patient demographics 


Characteristic 


IVIean (range) or N 


Age (years) 


57.7 (22-83) 


Body weiglit (l<g) 


79.76 (32.8-154.6) 


Heiglit (cm) 


170.73 (146.0-210.8) 


Gender 






39 


Female 


37 


Hematocrit (%) 


35.0 (11.8-44.7) 


Blood urea nitrogen (mg/dl) 


14.9 (5-63) 


Total bilirubin (mg/dl) 


0.47 (0.1-1.8) 


SGOT (U/l) 


35.6 (10-138) 


SGPT (U/l) 


27.5 (3-249) 


Alkaline phosphatase (U/l) 


178.7 (35-1,119) 


Creatinine (mg/dl) 


0.85 (0.3-1.4) 


GFR (ml/min/BSA) 


88.69 (42->120) 


Albumin (g/dl) 


3.80 (2.5-4.6) 


WBC (K/|jl) 


7.00 (2.1-27.8) 


Red blood cells (M/|jl) 


3.975 (2.27-7.18) 


Hemoglobin (g/dl) 


11.83 (7.5-15.7) 


Platelets (K/|jl) 


293.3 (75-896) 


Total protein (g/dl) 


6.91 (5.3-9.0) 


Triglycerides (mg/dl) 


113.4 (44-514) 


Glucose (mg/dl) 


105.9 (51-213) 


Cholesterol (mg/dl) 


175.7 (94-286) 



BSA, body surface area; GFR, glomerular filtration rate; SGOT, serum 
glutamic oxaloacetic transaminase; SGPT, serum glutamic pyruvic 
transaminase; WBC, White blood cells. 



Table 3 Treatment and sampling schedule for the clinical trials 









Dose 






Trial 




Number of 


range 


No. of 


PK sampling time 


No. 


Treatment 


patients 


(mg) 


observations 


points (days) (h) 


1 


Once 
weekly 


13 


10-60 


41 


0, 2,4, 24, 76, 168 
(day 1 , 2, 4, 8, 25, 
29, and 36)^ 


2 


Once 
weekly 


36 


1-16 


619 


0, 0.25, 0.5, 1,2, 
4, 8, 24, 48, 168 
(dayl) 


3 


Once 
weekly 


19 


15-35 


368 


0, 0.25, 0.5, 1,2, 
4, 8, 24, 48, 168 
(dayl) 


4 


Once 
daily 


8 


4 


98 


0, 1, 4, 6, 8, 24 
(day 14 and 28) 



PK, pharmacokinetic. 

^Sampling schedule is not same among the patients in trial 1 . 



the objective function value) for the fonA/ard addition and 0.01 
(equal to a decrease of 6.63 in the objective function value) for 
the backward elimination for 1 degree of freedom. 

Interindividual variability was described with a statistical 
model. Pharmacokinetic parameters were expressed by the 
exponential Eq. 1: 

Pij=PTVj-^M^ij) 

where P. represents the fth basic pharmacokinetic parameter 
for the /th individual, Pj^ \s the typical population value of the 
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yth parameter, and r|.. is realization from a normally distributed 
interindividual variable with a mean of 0 and an estimated vari- 
ance of co^ , which is the deviation of P. from Ptvj- Residual 
unexplained variability was evaluated using additive, propor- 
tional, and their combined error models. The NONMEM code 
is provided in the Supplementary Data online. 

Model evaluation. The final model was evaluated using a non- 
parametric bootstrap and visual predictive check.^^'^^ The boot- 
strap evaluation involved resampling the original data set 1 ,000 
times (sampling with replacement). The median values and 
the 2.5th and 97.5th percentiles of the parameters obtained by 
this analysis were compared with the ones obtained from the 
covariance step in NONMEM from the original data set. The 
visual predictive check was generated using 1 ,000 simulations 
from the final model, for all dose levels in our study, to assess 
the predictive performance and to verify if the performance is 
consistent among the dose levels. A graphical comparison was 
made between observed data and the model predicted median 
and 95% prediction interval over time. 
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Study Highlights 

WHAT ISTHE CURRENT KNOWLEDGE ON 
THETOPIC? 

Previous pharmacokinetic studies of sirolimus in trans- 
plant patients reported linear pharmacokinetics. 

WHAT QUESTION DIDTHIS STUDY ADDRESS? 

In this study, we explored and described the nonlinearity 
of sirolimus using concentration data from patients with 
advanced solid tumors treated with a wide dose range 
(1-60 mg/week). We developed a PopPK model with 
saturable absorption characteristics. 

WHATTHIS STUDY ADDSTO OUR KNOWLEDGE 

This is the first PopPK report indicating the nonlinearity 
of sirolimus and the first PopPK report in patients with 
cancer. 

HOWTHIS MIGHT CHANGE CLINICAL 
PHARMACOLOGY ANDTHERAPEUTICS 

Our model can be used to predict the effect of dose and/ 
or schedule change on the whole blood concentrations 
of sirolimus. This will facilitate future studies of this agent 
inaddition to providing information that may guide dosing 
for patients individually. 
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